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, Abstract. There are two fundamental problems studied by the theory of hamiltonian 

integrable systems: integration of equations of motion, and construction of action-angle 
variables. The third problem, however, should be added to the list: separation of vari- 
ables. Though much simpler than two others, it has important relations to the quantum 
' integrability. Separation of variables is constructed for the SL{3) magnetic chain — an 
example of integrable model associated to a nonhyperelliptic algebraic curve. 
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: 1 Introduction 

> . 
• 1—1 , 

^ ■ Consider a completely integrable Hamiltonian system with D degrees of freedom. 
. According to the definition of complete integrability due to Liouville-Arnold it 
means that the system possesses exactly D independent Hamiltonians Hj commuting 
with respect to the Poisson bracket 

{H^,Hk} = j,k = l,...,D (1) 

There are three fundamental problems discussed in the theory of integrable sys- 
tems. They are listed below in the order of decreasing complexity: 

• Construction of action-angle variables. 

• Integration of equations of motion. 

• Separation of variables. 
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For the wide class of finite-dimensional integrable systems subject to the Inverse 
Spectral Transform Method an effective integration of equations of motion can be 
performed using the techniques of algebraic geometry 0. As for the effective con- 
struction of the action-angle variables, it is more difficult problem 0, especially 
when the reality conditions are carefully taken into account 0]. The case of systems 
associated to hyperelliptic spectral curves is studied in detail [0, Q , there being only 
some preliminary results for the non- hyperelliptic case f^. 

In the present paper the third mentioned problem (separation of variables) is 
studied. To be precise, the separation of variables is understood here as construction 
of D pairs of canonical variables Xj, pj {j = 1, . . . , D) 

{xj,Xk} = {pj,Pk} = {pj,Xk} = Sjk (2) 

and D functions $j such that 

H,,H2,...,Hd) = j = 1, 2, . . . , D (3) 

where Hj are the Hamiltonians in involution. 

The above definition is a paraphrase of the usual definition of separability of 
variables in the Hamilton- Jacobi equation Note that the canonical transforma- 
tion from the original variables to {xj,pj) may not necessarily be a pure coordinate 
change, as in textbooks on classical mechanics, but can involve both coordinates 
and momenta. 

The problem in question, being the simplest of the three, is rather neglected 
in the literature on the subject, though, in our opinion, it deserves attention at 
least for two reasons. First, the variables {xj,pj) serve usually as a raw material for 
constructing action-angle variables and integrating equations of motion. Second, the 
problem is interesting for the theory of quantum integrability, since the construction 
of separated variables usually has direct counterparts in the quantum case 0]. 

The construction of the variables {xj,pj) is well known for the case of the hyper- 
elliptic spectral curve 0, ^, |], |^ though its relation to the separation of variables 
is not always stated manifestly. The coordinates Xj are defined as the zeroes of the 
corresponding Baker- Akhiezer function, and the canonically conjugated momenta 
Pj, or sometimes exppj, usually turn out to be eigenvalues of the corresponding 
L-operator taken at the values of the spectral parameter equal to xj. The functions 
$j are then simply the characteristic polynomials of the L-operator. 

In the present paper we study the problem for the nonhyperelliptic case. The 
5*17(3) classical magnetic chain is chosen as a sample toy-model. Having in mind 
subsequent application to the quantum integrability we make extensive use of the 
classical r- matrix formalism We consider complexified version of the model 



in order to avoid the additional complications of the real case. Our construction 
of separated variables is quite elementary and does not involve any sophisticated 
algebro-geometric techniques. 



2 



2 Description of the model 



The model we are going to describe is the nonhomogeneous classical SL{N) magnetic 
chain. It is in a sense generic for the models related to the S'L(A^)-invariant classical 



For N = 2 the model was introduced in IITTI The continuous 



r-matrix 

version of the model was studied earlier in 
(Gaudin model) see Section 5. The quantum version of the model is well studied by 
means of the Bethe ansatz method 



13, 14, 15 . For a degenerate case 



16, m 18 



The model in question is described in terms of the variables 



M; J2a=i ^aa = 0) subjcct to the Poisson brackets 



J cM c(") 1 _ ( ci-rn) r 



(4) 



which define the Kirillov-Kostant Poissonian structure on the direct product of M 
orbits of the coadjoint action of the Lie group SL{N) on sl{N)*, see e.g. [|10|- It 
is well known that the center of the Poisson algebra is generated by the eigenvalues 
of the matrices 5^™) 



det(M + ^(")l 



N 

m 

a=l 



N 



U 



m) 



(5) 



0=1 



We shall assume that /i™) are fixed numbers. The Poisson bracket (|^) is thus 
nondegenerate on the manifold @ having dimension 2D = MN{N — 1) for the case 
of generic orbit (all eigenvalues of S^"^'^ are distinct). In what follows we always 
assume that the orbit is generic. 

Let Z be an invertible N x N number matrix having distinct eigenvalues, let 
(^^ = 1, • • • , M) be some fixed numbers, and m be a complex parameter (spectral 
parameter). Consider the product (monodromy matrix) 



T{u) = Z{u - 6m + 5(^^)) ...{U-62 + 5(2)) (m - 5i + S^' 



(6) 



Proposition 1 Matrix elements ofT{u) have the following quadratic Poisson brack- 
ets ^ 

{Ta,p^{u),T^2l32iv)} = ^^^(^«2/3i(w)^ai/32(^^) - ^ai/32 W^a2/3i M ) (7) 

The proof (see |jTO|) is based on the fact that the factors {u — 5m + S^"^^) have 
the same Poisson brackets (|^ which reproduce themselves for the product T{u) 
(Lie-Poisson group structure). 

1 2 

Using the notation ||T^ T = T®id, T = id®T one can put the formula (|^) into 
a compact form 

{Tiu),Tiv)} = -^[V,Tiu)Tiv)] (8) 
where V is the permutation operator in ® C^. 
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Let the spectral invariants ty{u) of the matrix T{u) be defined as the elementary 
symmetric polynomials of its eigenvalues 



U{u)=tTf\T{u), iy = l,...,N 

For example, 

h{u) = trT(M), hiu) = ^{ti^T{u) - tTT\u)), . . . 

tN^u) = detT(n) = d{u). 

Note that the central functions are contained in the determinant d{u) = 
det T(m), see @. 

Proposition 2 The non-leading coefficients at powers ofu of the polynomials ty{u), 
u = 1, . . . , {N — 1), form a commutative, with respect to the Poisson bracket 
family of MN{N — l)/2 independent Hamiltonians. 

Proof. The polynomial ti,{u) having power uM in u contributes uM Hamil- 
tonians (its leading coefficient is a number), the total number of Hamiltonians is 
M(l + 2 + ■ ■ ■ + (A^ - 1)) = MN{N - l)/2. The commutativity of t^{u) 

is a direct consequence of the fundamental Poisson bracket (Q), see |jTO[|. The in- 



dependence of the integrals of motion is proven in ^ for a different model 



(Gaudin model, see Section 5) but the proof is valid also for our case. Note that 
the assumption made concerning nondegeneracy of the spectrum of the matrix Z is 
essential for the independence of t,y{u). 

By virtue of the proposition and since the number of Hamiltonians constructed 
D = MN{N — l)/2 equals exactly half dimension of the phase space the system is 
completely integrable. Now we can turn to the problem of constructing the separated 
variables. 

Conjecture 1 There exist functions A and B on GL{N) such that the following 
two assertions are true. First, A{T) is an algebraic function andB(T), respectively, 
is a polynomial of degree D = MN{N — l)/2 of the matrix elements T^^. Second, 
the variables Xj, Pj {j = 1, . . . , D) defined from the equations 

B{T{x,)) = 0, P, = A{T{xj)) (9) 

have the Poisson brackets 

{x„ Xk} = {Pj, Pk} = 0, {Pj, Xk} = Pj6jk (10) 

and, besides, are bound to the Hamiltonians ty{u) by the relations 

det(P,-T(x,)) = (11) 

The last relation means simply that Pj is an eigenvalue of the matrix T{u) when 
u = Xj. Putting Pj = exppj we observe that (|TTp fits the form since the spectral 
invariants of T{u) contain only the integrals of motion. 

In the present paper we prove the Conjecture |l] for the cases N = 2 and = 3. 
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3 SL(2) case 

Though the construction of the separation variables for = 2 is described in 
we reproduce it here in order to fix notation and to prepare the discussion of more 
difficult = 3 case. 

The system has M degrees of freedom. The spectral invariants of T{u) are 

t{u) = ti{u) = trT(M), d{u) = t2{u) = det T{u) 

the trace t{u) containing M integrals of motion. 
Define A and i3 as 

A{T) = Tn BiT) = Ti2 (12) 

and Xj, Pj respectively by the formulas {^). For the polynomial B{u) = B{T{u)) 
to have M zeroes it is necessary that its leading coefficient Z12 be nonzero. It 
can always be done by a similarity transform QT{u)Q~^ which affects neither basic 
Poisson brackets (P), nor Hamiltonians t{u)^ since the matrix Z, by assumption, has 
nondegenerate spectrum. 

Since the matrix T{u) becomes triangular at m = Xj the quantity Pj is an eigen- 
value of T{xj) and satisfies therefore the secular equation (0) which in the two- 
dimensional case takes the form 

P/ - t{xj)Pj + d{xj) = j = 1,...,M 

Note that the secular equation defines a hyperelliptic algebraic curve relating Pj and 

Xj. 

To prove the Conjecture ^ it remains to calculate the Poisson brackets of P's 
and x's. 

Theorem 1 The Poisson brackets for Pj and Xj are given by (|7^j. 

Proof. Let A{u) = A{T{u)) and B{u) = B{T{u)). Taking particular values of 
indices in (0) one obtains the identities 

{Aiu),Aiv)} = Q (13) 

{Biu),B{v)} = (14) 

{A(.),P(.)}^^i!M^l^^M^ (15) 

U — V 



The commutativity of P's ([1^ ) entrains obviously the commutativity of Xj (ze- 
roes of B{u)). The Poisson brackets including Pj can be calculated using implicit 
definition of Xj. From B{xj) = it follows that 

= {P,P(x,)} = {F,Biu)}u=., + B\xj){F,Xj} 

or 

{F,B{u)}.=,^ 

{F,Xj} = - 



B'{xj] 



for any function F. In the same way we have 

{P„F} = {A{x,),F} = {A{u),F}^=^^+A'{x,){x„F} 

Now it is easy to prove that {Pj,Xk} = Pj^jk- Starting with 

{Pj,Xk} = {A{u),Xk}u=x, + A'{xj){xj,Xk}, 

expanding the first term further, noting that the second term is aheady shown to 
vanish (0), and using (|I3p we arrive at 

{Aiu),Biv)}u=., ^ ^ B{x,)Aixk) - Aix,)Bixk) 
^ ""^"^ B'ixk) Xj-Xk B'{xk) 

The last expression vanishes for Xj 7^ Xk due to B{xj) = B{xk) = and is 
evaluated via L'Hopital rule for xj = Xk to produce the proclaimed result. The 
commutativity of P's can be shown in the same way starting from (|T3|). 



4 SL(3) case 

Let now N = 3. The polynomial T{u) takes values in 3 x 3 matrices. 



Tin) 



/ Tii(u) Ti2(m) Ti3(u) 
T2i(u) T22(m) T23(n) 
V T3i{u) T32(m) T33(n) 



The system has D = 3M degrees of freedom. The spectral invariants of the 
matrix T{u) 

ti{u) = trT(M) = A1 + A2 + A3 

h{u) = i(tr2T(M)-trT2(M)) = A1A2 + A1A3 + A2A3 

d{u) = det T('u) = A1A2A3 

are the coefficients of the characteristic polynomial for T{u) 

det(A - T{u)) = A^ - ti{u)X'^ + t2{u)X - d{u) 

which defines a nonhyperelliptic algebraic curve. 

It is convenient to introduce the matrix U{T) for any T G GL{3) 

U{T) = TAT = {T-ydetT 

/rri rri rri rri rri rri rri rri rri rri rri rri 

22-' 33 ~ -'23-' 32 -'23-'31 ~ -'21-'33 21-' 32 ~ 22-' 31 

13-^32 "~ -'12-' 33 -' ll-'33 ^ -' 13-'31 -' 12-' 31 ~ -' 11-' 32 

rri rri rri rri rri rri rji rji rri rri rji rji 

\ -' 12-' 23 ~ -' 13-' 22 -' 13-' 21 — -' 11-' 23 -' 11-' 22 " -' 12-' 21 

whose elements Ua/3 are algebraic adjuncts of T^,^. 



6 



Let U{u) = U{T{u)). The Poisson brackets for T and U are calculated easily 
from (1): 

{Tiu),U{v)} = — [V^ ,T{u)U{v)] (16) 



u — V 



{Ta,p,{u),Ua^f3^{v)} 



I](^^"i"2^7/3i(^)^7/32(^) + Ta^j{u)Ua2jiv) 5(3^(3^) (17) 

{U{u),U{v)} = [V, U {u)U{v)] (18) 

u — V 



(the superscript t2 in (|1^) denotes the transposition with respect to the second space 
in C^^C^^). 

The experience of the Inverse Spectral Transform Method and, in particular, 
SL{2) case suggests that in S'L(3) case the separated coordinates Xj, j = 1, • . . , 3M 
should be defined as zeroes of some polynomial B{u) of degree 3M and the corre- 
sponding momenta pj should be bound to Xj by the secular equation 

Pj - ti{xj)Pf + t2{xj)Pj - d{xj) =0, Pj = exppj 

It means that Pj should be an eigenvalue of the matrix T[xj). Therefore, there 
must exist such a similarity transformation 

T{xj) — > f{xj) = KjT{xj)Kr^ 

for each j that the matrix T{xj) is block-triangular 

fuixj) = fi3(a;,) = (20) 

and Pj is the eigenvalue of T{xj) splitted from the upper block 

Pj=f nix,) (21) 

The problem is reduced thus to determining the polynomial B{u) and the ma- 
trices Kj. Let us take the simplest possible triangular, one-parametric matrix K{k) 

K{k) = 

Note that the matrix 

f{u,k) = K{k)T{u)K-\k) 
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depends on two parameters: u and k. Therefore, we can consider the condition (^) 
as the set of two algebraic equations 

Ti2{x,k) = Ti2{x) + kT22{x) - kTii{x) - k'^T2i{x) = 

T,3{x,k) = T,3{x) + kT2s{x) =0 ^ ' 

for two variables x and k. Eliminating k from (^) one obtains the polynomial 
equation for x 

T23(Ti2T23 — T13T22) — Ti3(Ti3T2i — T11T23) = (23) 



or 



T23(x)f/3l(x) - Ti3(x)[/32(X) = (24) 

Since the matrix Z is assumed to have simple spectrum, the leading coefficient 
B{Z) of the polynomial B{T{u)) can always be made nonzero by a similarity trans- 
formation QT{u)Q~^, the equation (24) being thus of degree 3M. 



Expressing k from T13 = as A; = —Tis{x)/T23{x) and substituting it into the 
definition (|2T|) of P one obtains 

P = Tn{x) + kT2,{x) = (25) 

So, we have constructed 3M pairs of variables Xj, Pj. To prove the Conjecture 
1 it remains to show that they have good Poisson brackets. 

Theorem 2 The Poisson brackets for xj and Pj are given by (|7^). 

Proof. Let 

A{T) = JMH B{T) = T23Usi{T) - ri3W32(T) (26) 
-'23 

Putting A{u) = A{T{u)), B{u) = B{T{u)) and using (^, (0), (||) one can easily 
calculate the following Poisson brackets: 

{A{u),A{v)} = {B{u),B{v)} = (27) 



{A{u), B{v)} = (A{u)Biv) - B{u)A{v)^^] (28) 

from which the wanted Poisson brackets for xj and Pj are derived immediately in 
the same manner as in the SL{2) case. 

Remark. As N. Reshetikhin pointed out to us, the polynomial B(T), see (p6D, is 
invariant under similarity transform QTQ~^ acting on first and second row/column, 
Q G SL{2) C SL{3). The S'L(2)-invariance of B{T) entrains invariance of Xj and 
Pj. The meaning of this fact is still unclear. 



8 



5 Gaudin model 



The above construction of separated variables can be applied also to another in- 
tegrable system — Gaudin model — which was introduced first in the quantum 
variant |2^, see also ||2^, Its classical version turned out to be a useful example 



for developing a general group-theoretic approach to integrable systems |jT9|, ^0 

[m 
^13 



The model is formulated in terms of the same SL[N) variables S^J^^ as in Section 



2, see (H). Consider the matrix function 



M c(m) 
m= 



where {6m}m=i some fixed parameters and ^ is a traceless number matrix having 
N distinct eigenvalues. In contrast with T{u), see (^), the matrix T{u) has linear 
Poisson brackets. 

Proposition 3 

{Tiu),Tiv)} = ^[V,Tiu) + Tiv)] (29) 



The proof is a matter of direct computation [|l0 



Consider now the spectral invariants Ti,{u) of the matrix T{u) 

T^{u) =tiT{uY, iy=l,...,N 
Note that Ty{u) is a meromorphic function of u 



M u _a 



m=l a=l V"^ 



Note that C = trZ" and r^_^ = tr[S("^)]'^ = E^=i[4"'^]'' are numbers, see (g). The 
following proposition, analogous to Proposition 2, states the complete integrability 
of the system. 

Proposition 4 The quantities t!^^, (m = 1, . . . , M; u = 2, . . . , N ; a = 1, . . . , (z/ — 

1)), form a commutative, with respect to the Poisson bracket family of MN{N— 
l)/2 independent Hamiltonians. 

Proof. It is easy to compute the total number of Hamiltonians M(l + 2 + ■ ■ ■ + 
(A^ — 1)) = MN{N — l)/2 which equals exactly half dimension of the phase space. 
The commutativity of the spectral invariants of T{u) follows directly from (p9D, see 
T9| , The independence of the Hamiltonians is also proven there. 



The analog of Conjecture 1 for the Gaudin model is presented below. 

Conjecture 2 Let A and B he the same functions on GL{N) that in the Conjecture 
1. Then the variables Xj and pj defined by the equations 

B{T{x,)) = 0, p, = Ainx,)) (30) 

have the canonical Poisson brackets ^ and, besides, are bound to the Hamiltonians 
'^mu by relation det{pj — T{xj)) = 0. 
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The separation of variables for N = 2 and = 3 cases is performed now in the 
same manner as, respectively, in Sections 3 and 4. 



Theorem 3 The Conjecture 2 is true for N = 2. 

Proof. In the SL{2) case, the functions A and B on GL{2) are defined by 



the formulas ([T2|). Like in Section 3, we can always suppose that B{Z) ^ 0. The 
variables xj, pj, j = 1, . . . , N are then determined by the equations (pO|). 

Let Ag{u) = A{T{u)) and Bg{u) = B{T(u)). Taking particular matrix elements 



of ( P9[ ) one obtains 

{Ag{u), Ag{v)] = {Bg{u), Bg{v)] = (31) 

{Ag{u), Bg{v)} = (32) 

u — V 

The rest of the proof follows that of the Theorem 1. 
Theorem 4 The Conjecture 2 is true for N = 3. 



Proof. In the SL{3) case define A, B by the formulas (|26D and, like in Section 
4, suppose that B{Z) ^ 0. 

Let again Ag{u) = A(T{u)) and Bg{u) = B{T{u)). It suffices to estabhsh the 
Poisson brackets 

{Ag{u), Ag{v)} = {Bg{u), Bg{v)} = (33) 

{Ag{u\ Bg{v)] = ^ (bg{v) - Bg{u)^^^\ (34) 

since the remaining calculation is standard. 

It is possible to verify the above Poisson brackets directly, using (pQ]). It is 



simpler, however, to avoid long computations and to use the fact ||2l|, ^ that the 
Gaudin model is in fact a degenerate case of the magnetic chain. To be precise, let 
us replace S^"^^ in (D by Z by 1 + sZ, and divide T{u) by 11™=! (« - ^m)- 

Then, in the first order in e, we have T{u) / 11™=! {u — 5m) = 1 + sT{u) + O {e^) . The 
Poisson brackets ( p9D are obtained, respectively, as the linearization of the quadratic 
Poisson brackets (|). 

To conclude the proof, it remains to notice that 



A{u) = 1 + eAG{u) + 0{e') B{u) = e^BG{u) + 0{e*) 

and that the Poisson brackets (p3D, (p^) are obtained in the leading order in e from 
271), dH). 
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6 Unsolved problems 



The natural question arises whether the construction of separated variables pre- 
sented here for the SL{2) and S'L(3) cases can be generalized to the SL{N) case 
and further, to the integrable systems associated with classical r-matrices corre- 
sponding to other simple Lie algebras. Hopefully, the generalization will elucidate 
the geometric and algebraic meaning of the construction. 

The SL{N) case is presently under study. The problem consists in finding a 
multiparametric family of matrices K{ki, . . . , kg) such that after eliminating /c's 
from the system Ti2{x) = . . . = Ti]\f{x) = the resulting equation for x provide 
the necessary number of commuting zeroes. Another challenging object of study is 
the Kowalewski top which can be considered as a Gaudin model for 5^(4) ^ 50(5) 



group p4 |. 



Since the construction of separated variables for the SL{2) case has the direct 
quantum counterpart [|, 0, §], it seems reasonable to conjecture that the same is 
true for the SL{N) case. 

Acknowledgments. I am grateful to A. Bobenko, J. Hietarinta, V. Kuznetsov, N. 
Reshetikhin and A. Reyman for valuable and encourageing discussions. 
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